The Fergana Valley Is an Isolate of Biodiversity: A Discussion of the Endemic Herpetofauna and Description of Two New Species of Alsophylax (Sauria: Gekkonidae) from Eastern Uzbekistan

Simple Summary This study was carried out in one of the most densely populated and geographically isolated regions in the Republic of Uzbekistan, the Fergana Valley. The Fergana Valley has the highest level of endemic biodiversity in Uzbekistan (and one of the highest in Central Asia), and the habitats of these endemic species are rapidly being developed for agricultural purposes. Given this development, the remaining areas of habitat are not being adequately protected. The main goal of this study was to obtain up-to-date data on the distribution and abundance of five endemic reptile species in the remaining isolated and undeveloped habitats across the Fergana Valley. One of the most important achievements was the discovery of two unique and new micro-endemic species of gecko genus Alsophylax, which are described herein. These results elevate the number of endemic species in the Fergana Valley and further highlight the urgent need to create state-protected areas of habitat with IUCN I and II protection status for the remaining areas of suitable habitat, which is currently not available. Abstract The high level of endemism in Fergana Valley has been well documented in numerous studies for various groups of animals and plants. In a relatively small area, there are 45 endemic plant species, five endemic insect species, and five endemic reptile species. In surveying this area for data on distribution, abundance, acoustics, and genetic samples for species of reptiles, we discovered two new species of gecko from the genus Alsophylax. Phylogenetic analyses of mitochondrial DNA sequences indicate the relatives of these new species are the even-fingered gecko, Alsophylax pipiens, and the southern even-fingered gecko, Alsophylax laevis, located hundreds of kilometers to the northwest and southwest of the Fergana Valley. The threats to these new endemic species are significant given the amount of continued agricultural development that involves new territories previously considered “unsuitable” for any species of significance that is leading to the further reduction in, fragmentation of, and degradation of the remaining natural ecosystems in the Fergana Valley. The conservation of these rare and locally endemic species depends directly on the readiness of the state to create areas with IUCN I and II protection. The many studies documenting levels of endemism, along with the data published in this study, are the basis for the justification for state-protected areas in the Fergana Valley.


Introduction
The Fergana Valley is located in the eastern region of Uzbekistan and represents an ancient, isolated ecoregion with unique flora and fauna, much of which is endemic to the valley. In a relatively small area, the Fergana Valley contains several dozen endemic plant species, five endemic species of insects, and five endemic reptile species, most of which are associated with the sandy massifs and foothill habitats in the western and southern areas [1]. With the level of aggressive agricultural development and effects associated with climate change, each of these endemic species is facing its own set of environmental threats. For example, Strauch's toad agama (Phrynocephalus strauchi Nikolsky, 1899), Said-Aliev's toad-head agama (Phrynocephalus saidalievi Sattorov, 1981), Fergana sand racerunner (Eremias scripta pherganensis Szczerbak et Washetko, 1973), and Strauch's even-fingered gecko (Alsophylax loricatus Strauch, 1887) are all sandy massif specialists, and all five species fall into the "risk zone" [2]. These historically undisturbed sandy massifs, which comprise one of the key habitats for endemic species of reptiles in the Fergana Valley, have been reduced to 175.71 sq. km, and there are currently only 13 isolated plots of habitat, which now occupy an area of 0.46-58.66 sq. km [3].

The Current Status of Conservation in Fergana Valley
Agriculture and textile development has a long history in the Fergana Valley causing the gradual reduction in native habitats. The greatest reduction occurred in the early 20th century due to textiles and vegetable farming, leading to the development of the flatter areas such as sandy habitats. These sandy habitats represent the remains of a vast desert complex in the Fergana Valley originating from the Syrdarya River and are quite old, allowing for adaptation and speciation in these unique microhabitats. The rapid growth of the population of Uzbekistan [3] led to large portions of these habitats being plowed for development. Simultaneously, the rise in population also led to the construction of many agricultural fields, fish ponds, and houses, further fragmenting sandy habitats into sand islands, and preventing migration between adjacent populations. Additionally, excessive watering of adjacent agricultural developments has had a negative impact on existing habitat, resulting in a change to the hydrological regime [3]. This increase in groundwater creates more vegetation growth in these sand habitats, reducing the open cover dune habitat. This reduction in open dune habitat increases interspecific competition between the endemic sand-dune-adapted species and surrounding non-sand-dune-adapted species, making conservation efforts ineffective [4,5]. The problem of preserving sand massif habitats in the Fergana Valley is further complicated by the fact that establishing the correct status and protection regime for these unprotected areas in Uzbekistan is technically quite difficult.
Until 2023, conservation of these endemic reptile species and their habitats was hampered by the fact that, in the Fergana Valley, there were no pre-existing areas with IUCN I and II protection status. The existing natural state monument, "Yazyavan sands" [1,2], and other natural district-level monuments do not meet the requirement for higher-level biodiversity conservation [6]. However, since 2020, the project "Preservation of Key Natural Complexes in the Fergana Valley (Republic of Uzbekistan)" through the Michael Succow Foundation, with support from the World Wildlife Fund (WWF), Critical Ecosystem Partnership Fund (CEPF), and Central Asian Desert Initiative (CADI) project, has been garnering more attention for the conservation of the Fergana Valley. Additionally, the project "Maintenance and Conducting Selective Records of Rare and Endangered Species of Vertebrate Animals of the Fergana Valley", carried out by the "Cadaster and Cadastre of Rare Vertebrate Animals" at the Institute of Zoology of the Academy of Sciences of the Republic of Uzbekistan, has been documenting new populations of species in the Fergana Valley since 2020. In 2019, our team conducted a National Geographic Explorers expedition to survey reptiles in Uzbekistan including the Fergana Valley, and our team returned for multiple expeditions in the Fergana Valley between late 2019 and late 2022. Our team assessed the current state of the remaining preserved isolated habitat sites; as a result, we were able to update distribution and population density data for many reptile species in Uzbekistan. Most of our surveys were focused on the endemic and endangered species in Fergana Valley. Moreover, during these surveys, we discovered two new micro-endemic habitat specialist species of the small body-sized gekkonid genus, Alsophylax.

Current Taxonomy of the Gecko Genus Alsophylax
Alsophylax Fitzinger, 1843 originated from the description of Lacerta pipiens by Pallas [7] and was subsequently placed in its own genus, Alsophylax, by Fitzinger [8]. Since then, Alsophylax has received little attention due to the fact that its species are very small in size and have what can be perceived as a conserved morphology across the distribution of the genus. As a result, the taxonomic history of this group has been extremely complicated and confounding. Fortunately, there have been detailed descriptions of useful characteristics for describing new species in recent works [1,[9][10][11]. Currently, Alsophylax comprises six species distributed in Central Asia and southern Russia (A. laevis Nikolsky, 1907 [1]. To date, no study has investigated the interspecific relationships within Alsophylax, and even its phylogenetic position among other gecko genera remains poorly understood, further exemplifying their understudied nature [12][13][14].
In this study, we provide the most complete phylogeny for Alsophylax on the basis of the mitochondrial barcoding gene cytochrome oxidase I, and we include expanded sampling for species complexes such as A. laevis and A. pipiens. We used this phylogeny along with a dataset of morphological characteristics to demonstrate that these new populations of Alsophylax from the Fergana Valley are distinct species and clearly demonstrate the need for further work in Uzbekistan and the Fergana Valley to find other undescribed Alsophylax species. Lastly, we provide a discussion on the history of herpetological work in the Fergana Valley, highlighting its endemic reptile species and providing updated distribution and population density assessments of these endemic reptiles.

Population Density and Field Observations
We sampled individuals using the line transect survey method at each sand massif habitat island. During the field observations, we observed 1288 individuals of seven species along the total route length of 210 km. Diurnal species were estimated from routes with varying width along the transect. We recorded the perpendicular distance from the transect line to each individual reptile that was recorded to calculate the transect width We carried out nocturnal sampling using headlamps at a fixed width along the transect. When surveying Teratoscincus, we also used the counting of the reddish eye reflections with some specimens detectable at quite far distances up to 100 m or more. However, we understand that using a wide census band can underestimate occurrence data. Therefore, we limited it in accordance with the terrain and vegetation cover of the area. The accounting results, average detection distance (y) of the species and effective width of the census band (B), were calculated. We calculated population density (D) per hectare (ha) using the following formula: where N is the total number of individual reptiles, n is the number of individuals with measured detection distances, and L is the route length. Comparative analysis of the results obtained using various route accounting methods proved its high accuracy [15].
This method has been successfully employed in multiple studies to survey desert reptile species [16][17][18]. We calculated mean values of the population density along with estimation of their error (M ± m) according to the results at several approximately equal route sections. We conducted daytime surveys in sunny weather with low wind and no precipitation. We started recording data coinciding with the lizards' maximum activity, and finished when they shifted to shelters to escape the heat during the maximum daytime temperature. We also recorded substrate temperature.

Molecular Data and Phylogenetic Analyses
Our dataset included 46 samples belonging to all six Alsophylax species, including the two new species from the Fergana Valley. This also includes one Altiphylax species from Central Asia, which we hypothesize is an Alsophylax, and one sample of Mediodactylus kirmanensis (Nikolsky, 1900) from southern Iran as an outgroup. All analyzed materials are listed in Appendix A.
We extracted total genomic DNA from ethanol-preserved muscle or liver tissues using standard phenol-chloroform extraction procedures [19] followed by isopropanol precipitation. We amplified a fragment of the cytochrome oxidase I (COI) gene with maximal length of 655 bp. This mitochondrial marker is widely used for barcoding in vertebrates [20,21] and has also proven to be useful for species identification in reptiles [22][23][24]. Primers used both for PCR and sequencing were the VF1-d (5 -TTC TCA ACC AAC CAC AAR GAY ATY GG-3 ), VR1-d (5 -TAG ACT TCT GGG TGG CCR AAR AAY CA-3 ) [25], RepCOI-F (5 -TNT TMT CAA CNA ACC ACA AAG A-3 ), and RepCOI-R (5 -ACT TCT GGR TGK CCA AAR AAT CA-3 ) [24]. We sequenced fragments in both directions for each sample, and a consensus sequence was generated. We used 25 µL reactions for PCRs with approximately 50 ng of genomic DNA, 10 pmol of each primer, 15 nmol of each dNTP, 50 nmol of additional MgCl 2 , Taq PCR buffer (10 mM Tris-HCl, pH 8.3, 50 mM KCl, 1.1 mM MgCl 2 , and 0.01% gelatine), and 1 U of Taq DNA polymerase. We used the following PCR conditions: an initial denaturation step at 95 • C for 3 min; five cycles at 95 • C for 30 s, annealing at 45 • C for 1 min, and extension at 72 • C for 2 min; 35 cycles at 95 • C for 30 s, annealing at 48-51 • C for 1 min, extension at 72 • C for 2 min, final extension of 5 min at 72 • C, and storage at 4 • C. We loaded PCR products onto 1% agarose gels, stained with GelStar gel stain (Cambrex), and visualized in a Dark reader transilluminator (Clare Chemical). If results were satisfactory, products were purified using 2 µL, from a 1:4 dilution of ExoSapIt (Amersham), per 5 µL of PCR product prior to cycle sequencing. The 10 µL sequencing reaction included 2 µL of template, 2.5 µL of sequencing buffer, 0.8 µL of 10 pmol primer, 0.4 µL of BigDye Terminator version 3.1 Sequencing Standard (Applied Biosystems), and 4.2 µL of water. The sequencing reaction was 35 cycles of 10 s at 96 • C, 10 s at 50 • C, and 4 min at 60 • C. We purified cycle sequencing products using ethanol precipitation. We carried out sequence data collection and visualization on an ABI 3730xl automated sequencer (Applied Biosystems). Resulting sequences are deposited in GenBank under accession numbers (see Appendix A).
All Sequences were aligned manually using BioEdit Sequence Alignment Editor 5.0.9 [26]. The final alignment used for all phylogenetic analysis contained 654 bp of COI gene for all samples. MODELTEST v.3.06 [27] was used to estimate the optimal evolutionary models to be used for phylogenetic analysis. The preferred model was (GTR + G), as suggested by the Akaike information criterion (AIC). Confidence in tree topology was estimated with bootstraps in RAxML 8.1.1 [28]. We considered values above 85% as strong support and rather than choose an arbitrary number of bootstrap replicates to run we use internal metrics of the to allow the analysis to stop itself when sufficient replicates have been run.

Morphological Data and Statistical Analyses
We analyzed 310 specimens of all six recognized species of the genus Alsophylax for 22 morphological characteristics, including 14 metric and eight meristic characteristics. We chose these morphological characters on the basis of previous studies on Alsophylax [9], as well as additional new characteristics that were identified for this study. All specimens used were from Zoological Museum of Moscow State University (ZMMU), Zoological Institute of St. Petersburg (ZISP), Institute of Zoology National Academy of Science, Ukraine (IZ NAS), and Institute of Zoology Academy of Sciences of the Republic of Uzbekistan (IZ) (Appendix A).
Morphological measurements were taken with a digital caliper to the nearest 0.1 mm and included the following list of characteristics (14 in total): snout-vent length (SVL) measured from tip of snout to vent; tail length (TL), measured from vent to tail tip; head length (HL), distance between retroarticular process of jaw and snout tip; head width (HW), measured at the widest point of the head; head height (HeadH), measured from top of occiput to underside of jaws; orbital diameter (OrbD), diameter of orbit; snout-to-eye distance (SnEye), measured between anterior most point of eye and tip of snout; eye-to-ear distance (EyeEar), distance from anterior edge of ear opening to posterior corner of eye; humeral length (LS), measured on ventral surface of arm from base of axilla to posterior margin of elbow while forelimb bent by 90 • at elbow; forearm length (ForeaL), measured on dorsal surface while forelimb bent by 90 • at the elbow from the posterior margin of elbow to wrist; femur length (FemurL), measured on the ventral surface of the leg from base of femur to knee while hindlimb bent by 90 • ; crus length or tibia length (CrusL), measured on ventral surface of the of knee to the base of the heel while hindlimb bent by 90 • ; length of finger IV (LD4A), length of free distal phalanx of forth finger, without claw; length of toe IV (LD4P), length of free distal phalanx of forth toe, without claw.
The following characteristics (eight in total) were examined: number of scales around midbody including ventral and dorsal scales (Sq); number of scales along midbody from mental shield to anterior edge of cloaca (SLB); number of supralabials (SL); number of infralabials (IL); number of scales along middle of head, between eyes (I); number of subdigital lamellae under fourth finger (LF 4); number of subdigital lamellae under fourth toe (LT 4); number of postcloacal spurs on base of tail on both sides (Spurs); number of precloacal pores (PP).
We investigated morphometric differences between species of Alsophylax and the new species described herein. We log 10 -transformed all meristic measurements to normalize the data normality and increase the homogeneity of variance. We adjusted morphometric measurements to remove the effects of body size variation [29]. We use the allometric formula: Xadj = log(X) − b[log(BL) − log(BLmean)].We conducted a discriminant function analysis (DFA) to characterize the morphological variation between one new population of Alsophylax from the Fergana Valley and A. laevis and A. pipiens.

Population Density and Distribution
Through literary sources, as well as the results of our surveys of the Fergana Valley from 2019 to 2022, we noted 27 species of reptiles, which accounts for 43.5% of all reptile species known from Uzbekistan. Eight (30%) of the 27 Fergana Valley species are considered rare and listed in the Red Book of the Republic of Uzbekistan, five species (18.5%) are listed in the IUCN Red List, and three (11%) are CITES-protected species. Lastly, there are seven (26%) species that are micro-endemic species. The data from our surveys were analyzed and used to make Table 1, which lists all surveyed locations with an estimated average density of each endemic reptile species.

Molecular Analysis and Morphology
All phylogenetic analyses support the monophyly of Alsophylax ( Figure 1). These analyses demonstrate that both A. pipiens and A. laevis are polyphyletic, with significant genetic structure within each of these wide-ranging species complexes. These analyses also indicate that both of these new Alsophylax populations from the Fergana Valley ( Figure 2) are represent unique evolutionary lineages, sister to different clades of either A. pipiens or A. laevis. All sequences have been deposited on Genbank and their accession number can be found in Appendix B.
Taxonomy Family Gekkonidae Genus Alsophylax Fitzinger, 1843 Both new species belong to the genus Alsophylax according to the following morphological characteristics: (1) Elongate body with relatively short limbs; (2) non-segmented tail approximately the same length as the body; (3) relatively small and roundish head with a short rostral part; (4) no femoral pores and 6-8 precloacal pores; (5) small, roundish and flattened dorsal tubercles are present in some representatives of this genus, or dorsal scales can be smooth and homogeneous in another species.
Alsophylax ferganensis sp. nov. (Figures 3 and 4a).     Table  2).     Diagnosis. Alsophylax ferganensis sp. nov. is tentatively a sister to clades A of A. pipiens and A. laevis ( Figure 1) and morphologically closer to A. laevis (Table 3). Alsophylax ferganensis sp. nov. sp. nov. can be distinguished from A. laevis by a smaller maximum body size (SVL max 31.5 mm versus 38.7 mm) and relatively narrow head, as well as elongated limbs and dorsal patterns with well-defined nuchal loop, and relatively narrower dark transverse bands with approximately equal interspaces in between versus wide transverse patterns and narrow interspaces in A. laevis ( Figure 6). Caudal margins have dark transverse bands that are wavy. Dorsal scales are flat, smooth, and roundish, without enlarged dorsal tubercles (Figure 7).  Diagnosis. Alsophylax ferganensis sp. nov. is tentatively a sister to clades A of A. pipiens and A. laevis ( Figure 1) and morphologically closer to A. laevis (Table 3). Alsophylax ferganensis sp. nov. sp. nov. can be distinguished from A. laevis by a smaller maximum body size (SVLmax 31.5 mm versus 38.7 mm) and relatively narrow head, as well as elongated limbs and dorsal patterns with well-defined nuchal loop, and relatively narrower dark transverse bands with approximately equal interspaces in between versus wide transverse patterns and narrow interspaces in A. laevis ( Figure 6). Caudal margins have dark transverse bands that are wavy. Dorsal scales are flat, smooth, and roundish, without enlarged dorsal tubercles (Figure 7).    (Table 3).  Description of holotype. This is an adult male SVL 29.6 mm with a small teardrop shaped head (HL/SVL 0.27), smoothly passing into the neck. The rostral part is roundish and short (SnEye/HL 0.37), more elongated than the occiput (SnEye/EyeEar 1.11). The ear opening is small and roundish. There is a relatively large eye (Orb D/HL 0.25). The distance between posterior margin of the eye and ear (EyeEar/HL 0.33) exceeds the diameter of the orbit of the eye. Along the midline of the rostrum, a longitudinal depression is noticeable. The nostril is surrounded by the first supralabial, rostral, and enlarged nasal scale. Seven supralabials and five large infralabials of a rectangular shape Description of holotype. This is an adult male SVL 29.6 mm with a small teardrop shaped head (HL/SVL 0.27), smoothly passing into the neck. The rostral part is roundish and short (SnEye/HL 0.37), more elongated than the occiput (SnEye/EyeEar 1.11). The ear opening is small and roundish. There is a relatively large eye (Orb D/HL 0.25). The distance between posterior margin of the eye and ear (EyeEar/HL 0.33) exceeds the diameter of the orbit of the eye. Along the midline of the rostrum, a longitudinal depression is noticeable. The nostril is surrounded by the first supralabial, rostral, and enlarged nasal scale. Seven supralabials and five large infralabials of a rectangular shape are present, which decrease in length in the caudal direction. A large pentagonal mental shield borders on the first labial and two rounded postmentals. There are 18 round-oval scales between the eyes. The upper part of the head is covered by round, small granular scales which are enlarged in the area between the orbits. Large rounded gulars gradually pass into ventral large flat ventral scales (Figure 5b-e).
The body is elongated and slightly flattened, with no lateral folds. The dorsal scales are relatively small compared to the ventral; slightly larger and smaller scales are evenly distributed among them. No enlarged dorsal tubercles. A total of 86 scales are present from postmentals to cloaca, along with 44 around the midbody. Seven elongated, rectangular, perforated precloacal scales form an inverted V-shaped row. The opening of one pore occupies half of the total area of the precloacal scale.
The forelimbs are slender, and the humeral length slightly exceeds the size of the forearm (ForeaL/LS 0.9). The hindlimb is elongated, stronger in structure than the front. The thigh and lower leg are of the same length (CrusL/FemurL 0.1). The forelegs are dorsally covered with smooth juxtaposed scales, while convex granular scales exist on the ventral side of the upper arm, overlapping each other on the forearm. The dorsal surface of the hindlimb is covered by the scales similar to the dorsal ones, in contrast to enlarged scales on the ventral side, smooth and imbricated. The fingers are slender and long. From above, the metacarpus and metatarsus are covered with scales passing from the forearm and lower leg, respectively; those from below the palm and sole are covered with small swollen rounded granules, and the lower surface of the fingers is covered with transverse subdigital lamellae. The fourth digit of the hindlimb is the longest, with 14 slightly swollen subdigital lamellae.
The tail is 36.9 mm longer than the body (SVL/TL 0.8). The base of the tail is slightly swollen; on each side, there are two pairs of postcloacal spurs. The roundish subcaudal plates on the medial line of ventral side of the tail are 1.5-2 times larger than the surrounding scales.
Coloration. The dorsum of the body, head, and limbs are light beige. The pattern on the dorsum is formed by 5-7 transversely dark brown bands with uneven edges; the interspaces between bands is approximately equal to or slightly greater than the width of the bands. There are 7-10 transverse dark stripes on the tail, about the same number as the width of the gaps. Patterns on the head are not pronounced, whereby there is a narrow nuchal loop. The dorsal surface of the limbs is covered with transverse dark bands and irregularly shaped spots. The ventral surface of the body is white.
Distribution. Alsophylax ferganensis sp. nov. is only known from the type locality on the border with Kyrgyzstan, on the southern edge of the Fergana Valley ( Figure 2). The distribution of A. ferganensis sp. nov. is probably limited to this small mountain range, and a preliminary estimation of the available habitat indicates there is less than 50 km 2 . During field observations, another similar population was found on Pap Adyrs, which is located on the opposite side of the Fergana Valley. This isolated population is located 75 km from the type locality of A. ferganensis sp. nov., and its taxonomic status remains unclear.
Habitat and natural history. Clay-variegated canyons with sandstone outcrops represent the habitat of A. ferganensis sp. nov. (Figure 4b). The vegetation in these habitats consists of various shrubs and other plants. Most of the geckos were found on open hills with puffy cracked soil and at the bottom of the canyon along the sides.
Comparison. All representatives of the genus Alsophylax are divided into two groupsone dominant group has enlarged dorsal tubercles (A. loricatus, A. pipiens, A. przewalskii, A. szczerbaki) and the other group does not (A. laevis, A. tadjikiensis). Alsophylax ferganensis sp. nov. belongs to the second group and can be distinguished from A. loricatus, A. pipiens, A. przewalskii, and A. szczerbaki by granular homogeneous dorsal scales without enlarged tubercles ( Figure 6). A. ferganensis sp. nov. sp. nov. is morphologically closest to A. laevis but differs further by smaller maximal body size, relatively elongate limbs, and a relatively smaller head that is not sharply detached from the body; dorsal patterns are relatively narrower dark transverse bands with approximately equal interspaces in between, in contrast to the well-defined nuchal loop in A. laevis with wide transverse dorsal patterns and narrow interspaces. Alsophylax ferganensis sp. nov. differs from A. tadjikiensis by the following features: 6-8 precloacal pores in new species versus 5-7 pores in A. tadjikiensis; transversal dorsal bands versus pattern-less dorsum; bright yellowish or orange coloration of tail in A. tadjikiensis.
Etymology. This species is named after the Fergana Valley where it is endemic. Alsophylax emilia sp. nov. (Figures 3 and 8a) Animals 2021, 11, x FOR PEER REVIEW 17 of 28 Diagnosis. Alsophylax emilia sp. nov. is a sister to the clade containing lineages of A. pipiens, A. laevis, and A, ferganensis sp. nov. (Figure 1). Alsophylax emilia sp. nov. is   Diagnosis. Alsophylax emilia sp. nov. is a sister to the clade containing lineages of A. pipiens, A. laevis, and A, ferganensis sp. nov. (Figure 1). Alsophylax emilia sp. nov. is morphologically closest to the A. pipiens (Figure 3; Table 3) and can be distinguished by the following features: enlarged dorsal tubercles flat roundish smooth and larger than surrounding scales by no more than 2.5 times; head shape relatively more massive and sharply delimited from the body; 7-9 precloacal pores in V-shaped line; not elongated limbs; 5-6 wide transversal dorsal bands and narrow interspaces, no contrasting nuchal loop, and wide transverse patterns. The main distinguishing features of the new species are the following: maximum SVL of 35.0 mm; maximum TL 40.2 mm; 8/8 infralabials; 6/7 supralabials; two pairs of small roundish postmentals, which contact mentals; one nasal scale; 22 scales between the orbits; 97-104 longitudinal ventral scales from postmentals to cloaca; 59-65 scales along the midline around the body; 11-12 subdigital lamellae on the fourth finger; 13-14 subdigital lamellae on the fourth toe; 2-3 cloacal spurs on each side; medial line of subcaudal scales formed by noticeably enlarged plates, 2-3 times larger the surrounding scales (Table 4).
Description of holotype. This is an adult male, SVL 35.0 mm, TL 40.2 mm, with a relatively large head shape (HL/SVL 0.28) that tapers into the neck and body. The rostrum is rounded and relatively elongate (Snyer/HL 0.31). The distance between rostral and occipital parts is approximately the same (SnEye/EyeEar 1.03). The ear opening is small and oval. The eye is relatively small (Orb D/HL 0.22). The distance between the posterior margin of eye and anterior margin of ear (EyeEar 2.9) exceeds the diameter of the orbit (Orb D 2.15).
The nostril is surrounded by the first supralabial, rostral, and one enlarged nasal scale. There are eight supralabials and six large infralabials rectangular in shape, which decrease in size posteriorly. A large mental shield with roundish caudal edge borders with the first labial and four small and rounded postmentals. There are 22 granular scales between the eyes. The upper part of the head is covered by rounded scales which are enlarged on the sides of rostrum and between the orbits. Large, rounded gulars gradually pass into ventral large flat smooth juxtaposed scales (Figure 9).
The body is elongated and slightly flattened, with no lateral folds. The dorsal scales are heterogenic. Among small roundish scales, evenly distributed enlarged dorsal tubercles do not form regular rows (Figure 9d). Enlarged dorsal tubercles are no more than 2.5 times larger than surrounding scales. In the longitudinal medial row from postmentales to cloaca, there are 97 scales, along with 65 scales around the midbody. There are eight precloacal pores, located on elongated rectangular precloacal scales in an inverted V-shaped row ( Figure 9). The opening of the pore is large, occupying more than half of the total area of the precloacal scale. Between the row of precloacal pores and cloaca, there are five enlarged scales.
The forelimbs are slender, and the humeral length slightly exceeds the size of the forearm (ForeaL/LS 0.9). The hindlimbs are elongated, stronger in structure than the front. The thigh (6.5 mm) and lower leg are the same length. The dorsal surface of forelimbs and hindlimb is covered with smooth juxtaposed and imbricated scales, similar in size to the dorsal ones. The fingers are slender and not elongate (LD4A/SVL 0.07). From above, the metacarpus and metatarsus are covered with scales passing from the forearm and lower leg, respectively; those from below the palm and sole are covered with small swollen rounded granules, and the lower surface of the fingers is covered with transverse subdigital lamellae. The fourth digit of the hindlimb is the longest, with 13 slightly swollen subdigital lamellae.
The tail is 40.2 mm longer than the body (SVL/TL 0.9), cylindrical, rounded in crosssection, and not segmented. The base of the tail is slightly swollen; on the sides, there are two pairs of postcloacal spurs. The medial row features enlarged subcaudal plates.
Coloration. The background of the dorsum, head, and limbs are light beige (Figure 9d). The dorsal pattern is formed by six wide transversal dark brown bands with uneven edges; the interspaces between bands are approximately equal to or less than the width of the bands. There are 11 transverse dark bands on the tail. Patterns on the head are not pronounced, with a nuchal loop, dark brown strips on the sides of rostrum between eye, and supralabial plates. Moreover, no contrasting dark roundish spot is present on the tip of the rostrum. The dorsal surfaces of the limbs have transverse dark irregularly shaped spots (Figure 9). The color of the venter is white.
Distribution. Alsophylax emilia sp. nov. is known from a few localities on the northwestern edge of the Fergana Valley ( Figure 2). Specimens were found at the three localities on pap foothills (adyrs): vicinity of Jidalisai N40.891349; E70.884367; vicinity of Uygursay (N40.95126; E71.01833); 10 km northwest of Turakurgan village (N41.067757; E71.412176); one locality in the Namangan region, vicinity of Yartepa vill. N 41.067757, E 71.412176. We hypothesize that the distribution of A. emilia sp. nov. is probably broader in the Fergana Valley, but future studies are needed to clarify the distribution of this species in the surrounding regions.
Habitat and natural history. The habitat at all known localities of Alsophylax emilia sp. nov. was similar and made of clay canyons with outcrops, with sparse vegetation consisting of various shrubs (Figure 8b).
Etymology. This species is named in honor of Soviet and then Uzbek herpetologist Emilia V. Vashetko (24.04.1940-07.11.2022) for her great contributions to the study of the herpetofauna of Uzbekistan and surrounding countries.
Comparison. Alsophylax emilia sp. nov. belongs to the group of Alsophylax with heterogeneous dorsal scales and enlarged dorsal tubercles, which includes A. loricatus, A. pipiens, A. przewalskii, and A. szczerbaki (Table 3). Alsophylax emilia sp. nov. is morphologically closest to A. pipiens, from which it can be distinguished by the smaller dorsal tubercles on new species, which is are no more than 2.5 times smaller than surrounding scales; the medial line of subcaudal scales is formed by noticeably enlarged plates, which are 2-3 times larger than surrounding scales; there are 7-9 precloacal pores versus 8-13 in A. pipiens (Figures 8 and 9). Alsophylax emilia sp. nov. differs from A. loricatus, A. przewalskii, and A. szczerbaki by having smaller and smooth dorsal tubercles which are sporadically distributed across the back, versus large, strongly keeled trihedral dorsal tubercles organized in regular longitudinal rows; there are smaller scales on the head (22 scales between the eyes versus 8-15 in all three species mentioned above); dorsal patterns consist of transversal dark bands versus unclear dorsal patterns with longitudinal elements (Figure 6).
Alsophylax emilia sp. nov. further differs from A. tadjikiensis by the following characteristics: heterogeneous dorsal scalation versus homogeneous dorsal scalation; 7-9 precloacal pores in new species versus 5-7 pores in A. tadjikiensis; transversal dorsal bands versus pattern-less dorsum and bright yellowish or orange coloration of tail in A. tadjikiensis. Alsophylax emilia sp. nov. differs from A. ferganensis sp. nov. as follows: a larger body size with a maximum SVL of 35 mm versus 30.3 mm; heterogeneous dorsal scalation versus homogeneous dorsal scalation; elongate fore-and hindlimbs on A. ferganensis sp. nov.; transversal enlarged subcaudal plates versus roundish polygonal slightly enlarged subcaudal scales along the midline in A. ferganensis sp. nov. (Figure 10).

Taxonomy of Alsophylax and Conservation in the Fergana Valley
The genetic and morphological data in this study support the hypothesis that these new populations of Alsophylax represent undescribed new species, and that they currently appear to be endemic to the Fergana Valley of Uzbekistan (Figures 2 and 10). Given that

Taxonomy of Alsophylax and Conservation in the Fergana Valley
The genetic and morphological data in this study support the hypothesis that these new populations of Alsophylax represent undescribed new species, and that they currently appear to be endemic to the Fergana Valley of Uzbekistan (Figures 2 and 10). Given that these new species tentatively belong to separate species groups within Alsophylax suggests that the Fergana Valley may have been colonized by at least two different Alsophylax lineages (Figure 1). These results further highlight the age and complex geological processes that have helped generate the unique biodiversity in this isolated valley. Additionally, these new species elevate the number of endemic reptiles to the Fergana Valley and Uzbekistan. Lastly, we provided data on population densities for all endemic reptile species found in the Fergana Valley (Table 1; Figures S1-S4) indicating that not all these species respond to habitat degradation in the same manner. We discuss these results in conjunction with the history of previous herpetological studies in the Fergana Valley, and the future of conservation and species discovery in the area.

The History of Herpetological Study in Central Asia and the Fergana Valley
The history of studies on reptiles inhabiting the territory of Uzbekistan, as was the case with other vertebrates in Central Asia, began with the expedition of E. and gave a complete literary summary at that time [34]. The author processed the materials of A. P. Fedchenko and others who collected in Central Asia. The paper provides information on the distribution of seven species of amphibians, three species of turtles, 42 species of lizards, and 28 species of snakes. Of these, more than half were collected in Uzbekistan [34]. − Bogdanov published his work "Fauna of the Uzbek SSR. Amphibians and reptiles" [30] and provided all the data that he collected during his time doing field work in the Namangan region during 1949, 1950, 1954, and 1955.

Conclusions
Lastly, the age of the Fergana Valley allows it to support subtle yet different habitats ranging from claystone outcroppings to massive, isolated sand massifs and xeric flats that each support their own suite of species that are locally adapted to these different micro-habitats. The discovery of possibly two new species (Figures 1 and 3) that are microendemic habitat specialists isolated in various clay outcrops and canyons supports the hypothesis that the valley may have been colonized twice at different times. We hypothesize that this possible colonization pattern could indicate that other groups have experienced multiple colonization events and they too have levels of cryptic diversity among their populations. This study highlights the reptile diversity in the Fergana Valley and further supports the need for more surveys focusing on these habitats deemed "non-suitable", because the data from this study indicate that they may harbor unrecognized biodiversity. If the biodiversity from these areas is not documented and protected, they could become victim to the continued textile and agricultural developments in the Fergana Valley.

Acknowledgments:
The authors thank Ahmedov Askar, Soldatov Valentin, Adylova Dinara, Khasanov Mirpulat, Mun Yuriy, Bogatova Polina, and Chebotaryov Sergey for their assistance for collecting data in the field. They also thank the State Committee of the Republic of Uzbekistan for Ecology and Environmental Protection for their assistance in preparing permits for the collection of genetic samples.

Conflicts of Interest:
The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.
Appendix A Table A1. A list of specimens which were used in morphological and phylogenetic analyses.